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Abstract 

Here we develop a model equation of state which successfully parameterizes the thermody- 
namic functions of hadron resonance gas model at chemical freeze-out and which allows us to 
naturally explain the adiabatic chemical freeze-out criterion. The present model enables us to 
clearly demonstrate that at chemical freeze-out the resulting hadronic mass spectrum used in 
the hadron resonance gas model is not an exponential-like, but a power-like. We argue that 
such a property of hadronic mass spectrum at chemical freeze-out can be explained by the 
two new effects found here for wide resonances existing in a thermal environment: the near 
threshold thermal resonance enhancement and the near threshold resonance sharpening. The 
effect of resonance sharpening is studied for a sigma meson and our analysis shows that for the 
temperatures well below 92 MeV the effective width of sigma meson is about 50 to 70 MeV. 
Thus, the effect of resonance sharpening justifies the usage of the sigma-like field-theoretical 
models for the strongly interacting matter equation of state at such temperatures. Also we ar- 
gue that the most optimistic hope to find the quark gluon bags experimentally may be related 
to their sharpening and enhancement in a thermal medium. In this case the wide quark gluon 
bags may appear directly or in decays as narrow resonances that are absent in the tables of 
elementary particles and that have the width about 50-150 MeV and the mass about or above 
2.5 GeV. 

I. The hadron resonance gas model [1] is a reliable theoretical tool to extract information about 
the chemical freeze-out (FO) stage of the relativistic heavy ion collisions. However, the question 
about the reliable chemical FO criterion has a long history [H [21 13]. Very recently this question was 
thoroughly investigated again |3] , using the most sophisticated version of the hadron resonance gas 
model. Similarly to pj it was found that none of the formerly suggested chemical FO criteria [2] 
is robust, if the realistic particle table with the hadron masses up to 2.5 GeV is used. However, in 
[1] the criterion of the adiabatic chemical FO was suggested. In |1] it was also shown that despite 
an essential difference with the model used in [1] the same conclusion on the constant entropy 
per particle at chemical FO is well reproduced by the chemical FO parameters found in [H [5]. 
Furthermore, such a criterion of chemical FO is also consistent with the best description of the 
Strangeness Horn puzzle found recently in [6j . Thus, it turns out that the criterion of the constant 
entropy per particle at chemical FO is, indeed, the reliable one. It is interesting that the constant 
entropy per particle at chemical FO was also found in , but the way of hard core repulsion used in 
this work is too different from the traditional one used in the hadron resonance gas model [H SI E] 
and, hence, in contrast to the results of [H [3], the model used in leads to the simultaneous 
fulfillment of a few chemical FO criteria. 

Despite the long history of this question until recently there was no a single attempt to un- 
derstand what is the physical reason behind any of the chemical FO criterion. Only last year this 
problem got an adequate interest of theoreticians [3j. Here we develop an alternative approach and 
present a simple model equation of state which not only well describes the chemical FO thermody- 
namic parameters and reproduces the constant value of the entropy per particle at the chemical FO, 
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but in addition it allows us to elucidate the real mass spectrum of mesons and baryons that generate 
such a criterion. As it will be shown below the real mass spectrum of hadrons, i.e. the density of 
hadronic states, is very much different from the Hagedorn mass spectrum which is traditionally 
expected to emerge already for hadrons with masses above 1.2 GeV |H1 El HO]- However, below 
we demonstrate that the real mass spectrum of hadrons extracted from the adiabatic chemical FO 
model is not an exponential, but a power-like and the reason for such a behavior is an existence of 
many wide resonances. In this work we also show how at chemical FO a thermal medium essentially 
modifies the resonance mass distribution in case of large width leading to their narrowing near the 
threshold. Based on these findings we suggest that the quark-gluon bags may be observed at the 
NICA energy range as the narrow resonances (width about 50-150 MeV) with the mass about or 
above 2.5 which are absent in the tables of elementary particle properties. 

2. It is, indeed, a remarkable fact that at the chemical FO the entropy per particle deviates 
from 7.18 by about 8 % while the center of mass energy of collision changes from y^Sivw = 2.2 GeV 
to y/sjyN = 7 TeV in . In order to find the validity reason of the adiabatic chemical FO criterion we 
parameterize the mesonic (with subscript M) and baryonic (with subscript B) pressure as follows 



PM = Cm T^'' exp 
Pb = Cb T^^ exp 



Pm - rriM 



(1) 
(2) 



T 

>B - Mb 
. T 

assuming that the integrated mass spectrum of mesons and baryons can be represented by the 
constants Ca, Aa, ma, and ^m with a G {M,B}. Here T denotes the chemical FO temperature 
and fiB is the corresponding baryonic chemical potential. The pressure of antibaryons is, evidently, 
related to that one of baryons as Pb{T, Pb) = Pb{T, —Pb)- Note that Eqs. ([T]) and ^ represent the 
mixture of gases of massive particles with the temperature dependent number of degrees of freedom 
without the hard-core repulsion. The thorough inspection [U 0] shows that for the realistic values 
of hadronic hard-core radii below 0.45 fm the effect of hard-core repulsion can be safely neglected 
for temperatures below 150-170 MeV. This fact is implemented in Eqs. ([T]) and ([2]). 

Using the standard thermodynamic identities from Eqs. ([T]) and ^ one can find the particle 
densities pa = ^ = Pa/T with a G {M, B, B} and entropy density s = ^'^p^+vb+pb) chemical 
FO and get the following expression for entropy per particle 

s_ ^ (Am + ^^^™)pm + {Ab + M^){pb + Pb) - ^(Pb-Pb) 

Pp pM + Pb + Pb 

where the particle density is denoted as pp = pm + Pb + Pb- Iii order to qualitatively explain the 
reason of constant ratio s/pp we note that there exist two regions where the baryonic chemical 
potential roughly linearly decreases with the FO temperature IHH], i.e. Pb(T) ~ Pb(Tq) + p'b(T — 
To) ~ Pb{To) — p'b Tq + p'bT p^q + p'^ T. Taking the latter into account, one can write ([s]) as 



(Am + !Mm)p,, + {^Ab-p'b + '^^){pb + Pb) + 2^Pb 



(4) 



Pp Pm + Pb + Pb 

Consider, first, the chemical FO temperatures below 135 MeV [U H]. Note that in ^ the term 
<^ 1 since at large values of the baryonic chemical potentials the antibaryons are 
absent, while their number becomes comparable with the number of baryons for ^ — )■ 0. Now it is 
clear that one can get rid of the strong temperature dependence in (|4]), if the following conditions 
are obeyed 

rriM-pM Mb-P% 
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Then the entropy per particle at chemical FO can be cast as 



— ~ As - /i's - {Ab - ^i'B - Am) ^Ab- /i'fi , (6) 

Pp Pm + Pb + Pb 

i.e. up to a small correction the entropy per particle is constant. Strictly speaking, the above rough 
estimates should be valid only for low chemical FO temperatures, when the number of mesons is, 
indeed, a correction to the number of baryons. However, at high chemical FO temperatures the 
baryonic chemical potential is also linearly dependent on T |H H] , although with another value of 
the derivative jl'^. Then at high chemical FO temperatures the entropy density per particle (|5| 
becomes ^ ^ - /i'^ - {Ab - p'b - Am)^;^:^ ^ Am + {Ab - p'b - Am)^,^^^ ^ Const, 
since at these temperatures the fraction of baryons and antibaryons is almost a constant. Thus, the 
entropy per particle at chemical FO is basically defined by the power Aa of the dominated degrees 
of freedom {a = B for low T and a = M for high T) and by the speed of baryonic chemical potential 
change with the temperature change, i.e. by p'b- 

The above simplified consideration is well supported by a quantitative analysis. To demonstrate 
this in Fig. [l]we show the results of fit of the particle densities obtained in |4| for chemical FO. As 
one can see from Fig. [l] the quality of the chemical FO data for mesons (left panel) is much higher 
than that one of the baryons and antibaryons and, hence, the corresponding description of meson 
densities is excellent, while the description of baryon-antibaryon densities is good since the value 
of x^/dof = 7.8/11 is acceptable. The major part of this x^/dof value is generated by the three 
data points of the baryon-antibaryon density for the highest temperatures. We found the same 
minimum, if these three points are excluded from the fit, although in this case the corresponding 
value of x^/dof is much better /dof = 1.2/8. Our analysis shows that such a defect of the present 
model can be repaired by the excluded volume corrections. However, here we prefer to keep the 
model as simple as possible in order to demonstrate the main idea and, hence, hereafter we do not 
consider the last three data points into further fitting. 
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Figure 1: Temperature dependence of meson particle density pM (left panel) and baryon and an- 
tibaryon particle density (pb + Ps) (right panel) for the developed model. The data points with 
error bars are the results of the fit obtained in |3] for the experimental hadronic multiplicities. 



The results of the density fit are as follows 
Mb ~ 800.5 ±30 MeV, 
Cb ^ (2.624 ±0.191) ■ 10~^ 



MeV^^"^^ 



~ 6.097 ±0.38, A 



M 



fm^ ' 
5.31 ±0.14. 



-5 ± 5 MeV , 



Cm ^ (7.61 ±0.12) ■ 10" 



MeV 



fm'^ 



(7) 

(8) 
(9) 
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Note that the conditions ^ are well reproduced by the fitting of the densities. It may, however, 
look strange that for the vanishing value of mesonic chemical potential the meson mass parameter 
niM in ^ can be negative. We note that the results of fitting show that condition ([s]) for mesons 
is numerically well established, since the minimal value of the chemical FO temperature is about 
50 MeV. Therefore, the result ([T]) for mesons means that within the error bars one has hm = and 

UlM = 0. 

Using the above parameters we calculated (no additional fit!) the entropy density and found its 
ratio to the particle density. The results are shown in the left panel of Fig. |2} From this figure one 
can see that only the three points at highest chemical FO temperatures which were not included into 
a fit are not well reproduced. Of course, as it is seen from the right panel of Fig. [2] the additional 
fitting of the ratio s/pp can improve the quality description of all analyzed points and one gets 
xVdof ~ 0.47. 
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Figure 2: Left panel: Comparison of the chemical FO temperature dependence of the entropy per 
particle s/pp for the developed model (circles) and for the chemical FO data found in [3] (squares 
with error bars). Note that this ratio was just calculated using the parameters of Eqs. ([T]) - ([9]), 
i.e. no fitting was used. Right panel: the same ratio as in the left panel, but fitted together with 
particle densities. 

3. The model equations ([T]) and ^ allow us to directly find the mass spectrum of hadrons. The 
question of whether the experimental mass spectrum of hadrons given in the Particle Data Group 
tables coincides with the spectrum suggested by R. Hagedorn is of great interest nowadays [5|[U|[TU]. 
However, almost all discussions of the hadron mass spectrum simply ignore the width of resonances, 
whereas the large resonance width may essentially modify the spectrum [TTl [T2l [T3l ITll ITSl flGllTTllTS] 
and it may be even responsible for an absence of heavy excited resonances in the empirical mass 
spectrum of hadrons [HI |T9j. Therefore, it is interesting to study the effective mass spectrum of 
hadrons having the real width. Since in the hadron resonance gas model [U H] the resonance width 
is considered explicitly, its effect is also implicitly accounted for in the model equations Q and 

In order to elucidate the effective mass spectrum of baryons we rewrite their density in terms 
of the mass integral of the momentum integral of the Boltzmann distribution function. Since 
the effective baryon mass Mb in ([T]) is large compared to the maximal value of the chemical FO 
temperature, one can safely use the non-relativistic approximation for the momentum integration. 
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This means that the factor Ta should be assigned to the momentum integration, while the remaining 
T-dependence of the baryon (and antibaryon) particle density can be identically cast as 



PB 



exp 



Pb - Mb 
T 



C 



B exp 



Pb 
T 



dm 



Mb 



(m - Mb)^s 
T{Ab - 2.5) 



3.5 



T2 exp 



m 
T 
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{2tt)2Cb exp 



T 



dm 



(m - Mb) 



-3.5 



Mb 



ma 



T{Ab- 2.5) J (27r 



exp 



TFT 



(10) 



the momentum 



where in the last step of derivation we accounted for the fact that besides exp — ^ 
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integration in non-relativistic case generates the factor {2tt mT) 2 . Here T{A) is the usual gamma 
function and, hence, one has r(3.6) ~ 3.717. 



Since the original mass integral in (10) is the Laplace integral, then the representation (10) 



uniquely defines the T-dependence of baryonic density and vice versa. Thus, from ( 10 ) one concludes 



that at chemical FO the effective mass spectrum of baryons and antibaryons is rather a power-like 
than an exponential-like: 



dQBim) 
dm 



(27r)5 Cb {m-MB) 



r(3.6) 
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Figure 3: Ratio of the baryonic mass spectrum QbI^) of the present model found at chemical FO 
defined by (11) to the mass spectrum of all hadrons ^^^(m) taken from the Particle Data Group 
[20] and parameterized according to (12) [2T]. The full curve shows the ratio of integrated spectra 
QB{fn) / Q^B^{fn)i while the dotted curve demonstrates the ratio of the corresponding densities of 
states. The mass Mb = 800.5 MeV is taken from Eq. 

4. Comparing the baryonic mass spectrum f)B(m) = / dz -^M^ obtained from the density of 



states (11) with the parameterization 

exp / \ 

Qb {f^) 



m 



537 MeV 



Mr. 



5.72 



(12) 



suggested in [21] to describe the experimental hadronic mass spectrum, one finds a great disagree- 
ment between them! Indeed, as one can see from Fig. |3] these two mass spectra may deviate from 
each other by about 25 times at m ~ 1650 MeV and their asymptotic behaviors are completely 
different, since ^>s(m) ~ m^'^ in (11), while g^B^{m) ~ m^'''^ in (12). Naturally, there arise two 
questions, "Does it mean that one of these two mass spectra is wrong?' 



and "What is the reason 
for so huge difference between these mass spectra?" In order to demonstrate that none of these two 
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mass spectra is wrong and that at chemical FO the large difference between these two hadronic mass 
spectra is due to the width of hadronic resonances, we consider the Gaussian mass attenuation in- 
stead of the Breit-Wigner one that is used in the actual simulations since in this case the evaluation 
is more transparent. Also such a treatment would allow us to obtain some important conclusions 
on the mass spectrum of quark-gluon (QG) bags which according to [TSl HH] should unavoidably 
have the Gaussian mass attenuation. Note also that the estimates below provide us with the lower 
limit, since the Gaussian mass distribution vanishes much faster than the Breit-Wigner one. The 
typical term of the fc-resonance that enters into the mass spectrum of the baryonic particle density 
of the hadron resonance gas model is given by Fk{(Jk) exp ^ IH 15, where 



9k 



exp 



T 



(13) 



Here rrik is the mean mass of the fc-th resonance, Qk is its degeneracy factor, ak is the Gaussian 
width which defines the true width of such a resonance as Tk = Qak (with Q = 2 a/2 In 2) and the 
normalization factor is defined via the threshold mass Mj^ of the dominant channel as Nk{M'[^) = 



J dm exp 



2 at 



. For the narrow resonances the term Fk{(Jk) converts into the usual thermal 



density of particles, i.e. for 0"^ — )■ one has Fk — )■ gk 0(^fe, T), where the following notation is used 
0(m,T)^/(^expL^ 



A usage of Eq. (13) was heavily criticized in [21], but we find such a critique absolutely inad- 
equate for the states below the chemical FO. First of all, we note that in the approach of [23] the 
effect of medium cannot be switched off at any finite particle density or temperature. This means 
that according to the treatment of [23] (and many similar works!) all the hadrons whose momen- 
tum spectra are frozen due to the absence of any strong interaction between them should keep their 
momentum dependent width and mass which they acquired at the moment of kinetic FO up to they 
are captured by the detectors. Hence, according to [23] all hadrons measured by detectors, including 
the stable ones, are some resonances that 'feel' a thermal medium in which they were produced long 
after the medium is gone. Second, all the 'effects' which the authors of [24] claim to be of principal 
physical importance are reduced to a slight (by about 20 MeV) shift of the mass attenuation peak 
and small change of its shape for the A33 resonance, although the corresponding modification of the 
nucleon and pion properties the authors of [23] simply ignore. In our opinion such modifications 



of the mass attenuation of the A33 resonance compared to the 'crude' approximation of Eq. (13) 
cannot be measured in heavy ion experiments even for such narrow resonances as A33. Therefore, a 
serious discussion of similar 'effects' for heavy hadronic resonances whose mass and width are often 
known with the accuracy of 100 MeV or 200 MeV (or worse) [20] does not make any sense. Thus, 



at the present state of art there is no alternative to a physically transparent Eq. (13) to be used at 
and after the moment of chemical FO. 



2T7h^ 



The momentum integral in ( 13 ) can be written using the non-relativistic approximation 0(m, T) 

mT 



exp 



Then to simplify the mass integration of ( 13 ) one can make the full square in it 



rom the powers of (mfc — m) and get 
Fk{ak) = gk I dmfk{m) c:^ gk j dm j^^^^^Th^^ ' exp 



[ruk - m) 
2^^ 



mT 



exp 



mfc 
T 



,(14) 





where the following notations for an effective resonance degeneracy gk and for an effective resonance 
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mass ruk 



9k 



Qk exp 



2T2 



Qk exp 



k 



2Q2 T2 



mk = rrik- 



T 



^ fc 
Q2T 



(15) 
(16) 



are used. From Eq. ( 14 ) one can see that the presence of the width, firstly, may strongly modify the 



degeneracy factor gj. and, secondly, it may essentially shift the maximum of the mass attenuation 
towards the threshold or even below it. There are two corresponding effects which we named as the 
near threshold thermal resonance enhancement and the near threshold resonance sharpening. These 
effects formally appear due to the same reason as the famous Gamow window for the thermonuclear 
reactions in stars [^ : just above the resonance decay threshold the integrand in (14) is a 



product of two functions of a virtual resonance mass m, namely, the Gaussian attenuation is an 
increasing function of m, while the Boltzmann exponent strongly decreases above the threshold. 
The resulting attenuation of their product has a maximum, whose shape, in contrast to the usual 
Gamow window, may be extremely asymmetric due to the presence of the threshold. Indeed, as 
one can see from Fig. |4] the resulting mass attenuation of a resonance may acquire the form of the 
sharp and narrow peak that is closely resembling an icy slide. Below we discuss these two effects in 



some details. Qualitatively the same effects appear, if in ( 14 ) one uses the Breit- Wigner resonance 
mass attenuation instead of the Gaussian one. 



From the definitions of the effective resonance mass ( 16 ) and the effective resonance degeneracy 



(15) one can see that the effects of their change are strong for T ^ cr^. This can be clearly seen from 



Fig. |4| which demonstrates both of the above effects at low temperatures for the cr-meson. A simple 
analysis shows that the effect of resonance sharpening is strongest, if the threshold mass is shifted 
to the convex part of the Gaussian distribution in (14), i.e. for > fhk or for the temperatures T 



below 



mk-MTh — 13^- To demonstrate the effect of the width sharpening we list a few typical 
examples for baryons in the Table 1. For T < and for m > M'^^ the Gaussian mass distribution 



in ( 14 ) can be safely approximated as exp 



2-T 



exp 



Now recalling the standard definition of the width for the function /(x) = Q{x) Const exp [—hx 
one obtains the temperature dependent resonance effective width near the threshold as 



rf(T) 



ln(2) 

1 _ A 



ln(2) 



i_ 

T 



(17) 



since for such a distribution function /(x) one gets /(ln(2)/6) = /(0)/2. Note that in evaluating 
(17) we neglected the additional m^-^-dependence in (14), but one can readily check that numerically 
such a correction is small. The rightmost column in Table 1 demonstrates that Eq. (17), indeed. 



provides an accurate estimates for T < . The results of Table 1 also justify the usage of cx-meson 
and the field-theoretical models based on the well known cx-model for temperatures well below 

~ 92 MeV. Of course, the present approach which is developed for the chemical FO stage, when 
the inelastic reactions except for resonance decays are ceased to exist, cannot be applied for earlier 
stages of heavy ion collisions. However, here we would like to stress that an inclusion of the large 
width of (T-meson in the field-theoretical models of the strongly interacting matter equation of state 
is very necessary. From the above analysis one can see that the large width inclusion can generate 
some new important physical effects like the wide resonance sharpening in a thermal medium. 

From Fig. [sjone can see that the resonance enhancement can be, indeed, huge for wide (F > 450 
MeV) and medium wide (F ~ 300 — 400 MeV) resonances. This effect naturally explains the strong 
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Figure 4: Temperature dependence of the mass distribution fkijn) / (pim^^T) (in units of 1/MeV, 
see Eq. (14)) for cr-meson with the mass = 484 MeV, the width Fq- = 510 MeV [23] and 
MJ'^ = 2m7r — 280 MeV. In the left panel the short dashed curves below the two pion threshold 
(vertical line at m — mg- = —204 MeV) show the mass attenuation which does not contribute into 



the particle density (14). From the right panel one can see the effect of wide resonance sharpening 
near the threshold, i.e. an appearance of a narrow peak in the resulting mass distribution on the 
right hand side of the threshold which resembles an icy slide. For different temperatures this mass 
attenuation is shown by the solid, the short dashed and the long dashed curves. The a-meson 
effective width was found numerically from these mass attenuations: F^^^(T = 50 MeV) ~ 62.5 
MeV, VjfiT = 55 MeV) ~ 71.5 MeV and Vj^iT = 60 MeV) ~ 82.5 MeV. 



temperature dependence of hadronic pressure ([T]) at chemical FO, which in its turn generates the 
power-like mass spectrum of baryons (11). Clearly, the same is true for the mesonic pressure (|2| 
and mesonic mass spectrum in a thermal environment. However, we believe that a detailed study 
of such a phenomenon requires a special investigation. 

5. The first important conclusion from the analysis above is that there is no sense to discuss the 
mass spectrum of hadronic resonances, empirical or Hagedorn, without a treatment of their width. 
Furthermore, the same is true for the QG bags which, according to the finite width model [T9] . 
are heavy and wide resonances with mass Mb larger than Mq ^2.5 GeV and with the mean width 

of the form F^ ^ ro(T) ^ , where Fo(T) is a monotonically increasing function of T and Tq{T = 
0) e [400; 600] MeV. This range of Tq{T = 0) values corresponds to the cross-over temperatures 
Tco ~ 170 - 200 MeV [iHl HE] for vanishing baryonic density. The value Fo(T = 0) = 400 MeV 
is well consistent with the results of the present days lattice QCD thermodynamics [26|, [27j, but 
there is no guaranty that the lattice QCD data will not change in the future. Therefore, below we 
consider the whole range of values for the width Fo(T) analyzed in [TU] . 

There are two interesting features of QG bags which are related to the above treatment. Thus, 
from the results of [IHl HH] ^ind from (17) one can find the temperature Tg for the QG bags as 



T, 



B 



rg(T) 



'A/o(l-?fl) - (l-Ss) 



11.5-26. MeV, if 
46 - 104 MeV , if 
140 - 315 MeV , if 



ro 
ro 
ro 



0.4 - 0.6 GeV at T 
0.8 - 1.2 GeV at T 
1.4-2.1 GeV at T 



0, 

90, 

170. 



where = -jf^- denotes the ratio of the leading threshold mass M^^ of the bag to its mean mass 
Mb- In (18) the temperature T is given in MeV. Clearly for different bags the range of ^b value can 
be between and 1. Therefore, according to above results the bags with — > 1 should have been 



essentially enhanced and sharpened as the ordinary resonances. Moreover, according to (17) in this 
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Decay 
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62.5 
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-P33 
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A ^ ttA^ 


1080 


2.98 
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43.5 


N/A 


Pu 


1440 


350 


N ^-kN 


1080 


2.42 


38.74 


129.5 


N/A 


P33 


1600 


350 


A ttA 


1372 


1.53 
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68.7 


80.8 


P33 


1600 


350 


A ttA^ 


1080 
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30.3 


280. 


N/A 


Gn 


2190 


500 


pN 


1710 


2.26 


57.8 


74.6 


81.8 



Table 1: The parameters of several hadronic resonances together with their decay channels that 
are used to determine the quantities (3^ and . The last two columns show the corresponding 
effective width at temperature T = 50 MeV found, respectively, numerically from Eq. (13) and 



analytically from Eq. (17), when it can be applied. 
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Figure 5: Temperature dependence of the resonance enhancement. The ratio R{T) = ^^^^ 
is shown for the hadronic resonance decays given in Table I. For wide resonances the effect of 
enhancement can be huge. 



case for T ^ the QG bags should have had a small width F;^ ~ T+--r ^^(2) and, hence, such 
QG bags should have been stable or, in other words, these bags should have been observed! The 
reason why such bags are not observed in the experiments is naturally explained by the finite width 
model [IHllin]: it is due to the effect called as the subthreshold suppression, i.e. a huge suppression 
(of about fifteen to sixteen orders of magnitude compared to light hadrons!) of the QG bags for 
temperatures below the half of the traditional Hagedorn temperature Th (for more details see a 
discussion after Eq. (42) in [IH])- Such a suppression is a manifestation of the color confinement in 
terms of the QG bag width [T^ . 

On the other hand Eq. (18) also shows that the only hope to observe the QG bags exists, if 
1- Then for chemical FO temperatures much below Tg such bags could have sufficiently long 

eigen lifetime of about tb ^ ^ — 1^(2) — t in(2) • Substituting T ~ 0.5 T// G [80; 90] MeV in the 
last inequality and using the estimate of Eq. (18) for T = 90 MeV with ^b = 0.9, one finds the most 
optimistic estimates for the QG bag eigen lifetime as tb < 3.3 ± 0.3 fm/c. These estimates allow us 
to make the second important conclusion that an appearance of sharp resonances (baryonic or/and 
mesonic) with the width in the interval between 50 to 70 MeV aX the chemical FO temperatures 
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close to Tqgb — 85 ± 5 MeV that have the mass above 2.5 GeV and that are absent in the tables of 
particle properties would be a clear signal of the QG bag formation. Their possible appearance at 
chemical FO as metastable states of finite systems created in relativistic nuclear collisions is justified 
by the finite width model [ISl [IS]- At higher temperatures such QG bags can be formed too, but 
their width is larger and lifetime is shorter. The limiting values of at T for which the effect 
of resonance sharpening can exist is determined by the relation Tf^{T)/T > Q"^ Mo{l — C,b) from 
which one can see that the condition — )• 0.9 can be relaxed, but in this case the temperature of 
chemical FO gets higher. 

In addition one has to account for the statistical probability of the QG bags appearance at a 
given temperature T. Relatively to the nucleon the statistical probability of the QG bag of mass Mb 
is about W = ^'^ exp JIMn^Ms)] Rs{T), where Mn ^ 940 MeV is the nucleon mass and i?ij(r) 
is the resonance enhancement factor in a thermal medium. For T = 140 MeV and Mb = Mq ~ 2.5 
GeV one gets Wb — 3.85 ■ 1Q~^Rb- Our analysis shows that for such temperatures the typical 
resonance effective width values are about ~ 100 — 150 MeV while the typical values of the 
resonance enhancement factor is about Rb ~ 10 — 100 and, hence, compared to nucleon the relative 
statistical probability of such QG bags is about Wb — 3.85 ■ (10^^ — 10~^), which is essentially 
larger than the relative probability of the J /if) meson Wj/^ ^ 1.19 ■ 10"^ at the same temperature. 
Note that the chemical FO temperature T ~ 140 MeV corresponds to the highest SPS energy of 
collision at which the J /if) mesons are safely measured. Note that the chemical FO temperatures 
about Tqqb ~ 80 — 140 MeV correspond to the center of mass energy of collision ^Js^n ^ [4; 8] 
GeV [H H], which is in the range of the Dubna Nuclotron and NICA energies of collision. This 
energy range sets the most promising kinematic limit for the QG bag searches. 

6. In conclusion, here we present a simple model equation of state which successfully parame- 
terizes the thermodynamic functions of hadron resonance gas model at chemical FO. Such a model 
with the temperature dependent number of hadronic degrees of freedom allowed us to naturally ex- 
plain the adiabatic chemical freeze-out criterion that was found previously. In addition, this model 
allowed us to find out that the effective mass spectrum of baryons used in the resonance hadron 
gas model is not an exponential-like, but a power-like. Evidently, the same conclusion is valid for 
mesons, but such an analysis will be published elsewhere. 

In order to give a reason for the obtained difference between the effective baryonic mass spectrum 
and the empirical hadronic mass spectrum we analyzed the behavior of wide resonances in a thermal 
environment and found two new effects occurring, if the chemical FO temperature is small compared 
to the resonance width: the near threshold thermal resonance enhancement and the near threshold 
resonance sharpening. Further analysis showed that for the temperatures well below 92 MeV the a- 
meson can be rather narrow and it has the effective width of about 50 to 70 MeV. Thus, accounting 
for the (T-meson large width in a thermal medium allows us to justify the usage of the a-like field- 
theoretical models for the strongly interacting matter equation of state for temperatures well below 
92 MeV. Finally, we argued that the most optimistic hope to find the QG bags experimentally 
would be related to their sharpening and enhancement by a thermal medium. Then the QG bags 
may appear directly or in decays as narrow resonances of the width about 50-150 MeV that have 
the mass about or above 2.5 GeV and that are absent in the tables of elementary particles. 
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